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EXPLICIT FORMULAS FOR RELAXED DISARRANGEMENT DENSITIES 
ARISING FROM STRUCTURED DEFORMATIONS 
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Abstract. In this paper we derive explicit formulas for disarrangement densities of 
submacroscopic separations, switches, and interpenetrations in the context of first- 
order structured deformations. Our derivation employs relaxation within one mathe¬ 
matical setting for structured deformations of a specific, purely interfacial density, and 
the formula we obtain agrees with one obtained earlier in a different setting for struc¬ 
tured deformations. Coincidentally, our derivation provides an alternative method for 
obtaining the earlier result, and we establish new explicit formulas for other measures 
of disarrangements that are significant in applications. 
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1. Introduction 

Structured deformations provide a multiscale geometry that captures the contri¬ 
butions at the macrolevel of both smooth geometrical changes and non-smooth geo¬ 
metrical changes (disarrangements) at submacroscopic levels. For each (first-order) 
structured deformation {g, G) of a continuous body, the tensor field G is known to 
be a measure of deformations without disarrangements, and M :=Vg — G is known 
to be a measure of deformations due to disarrangements. The tensor fields G and 
M together deliver not only standard notions of plastic deformation, but M and its 
curl deliver the Burgers vector field associated with closed curves in the body and 
the dislocation density field used in describing geometrical changes in bodies with 
defects. Recently, Owen and Paroni I1T3I evaluated explicitly some relaxed energy 
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densities arising in Choksi and Fonseca’s energetics of structured deformations [^| 
and thereby showed: (1) {trM)'^, the positive part of trM, is a volume density of 
disarrangements due to submacroscopic separations, (2) {trMY , the negative part 
of trM , is a volume density of disarrangements due to submacroscopic switches and 
interpenetrations, and (3) \trM\, the absolute value of trM, is a volume density of 
all three of these non-tangential disarrangements: separations, switches, and in¬ 
terpenetrations. The main contribution of the present research is to show that a 
different approach to the energetics of structured deformations, that due to Baia, 
Matias, and Santos [1], confirms the roles of {trM)~^, [trM)~ , and \trM\ established 
by Owen and Paroni. In doing so, we give an alternative, shorter proof of Owen and 
Paroni’s results, and we establish additional explicit formulas for other measures of 
disarrangements. 

In order to motivate our study and to provide necessary background, we discuss 
briefly in the following subsections of this introduction some concepts and results 
from the multiscale geometry of structured deformations. (Readers familiar with 
this material may wish to skip to the last subsection of the introduction where our 
main results are summarized.) 

1.1. Structured deformations and disarrangement densities in the setting 
of Del Piero and Owen. The need in continuum mechanics to include the effects of 
multiscale geometrical changes led Del Piero and Owen 0 to a notion of structured 
deformations as triples (k, g, G) , where 

• the injective, piecewise continuously differentiable field g maps the points of 
a continuous body into physical space and describes macroscopic changes in 
the geometry of the body, 

• the piecewise continous tensor field G maps the body into the space of linear 
mappings on the translation space of physical space and satisfies the "accom¬ 
modation inequality" 

0 < C < det G{x) ^ det Vg{x) at each point x (1.1) 

where V denotes the classical gradient operator, and 

• K is a surface-like subset of the body that describes preexisting, unopened 
macroscopic cracks. 

A geometrical interpretation of the field G is provided by the Approximation The¬ 
orem DZI: for each structured deformation (k, g, G) there exists a sequence of injec¬ 
tive, piecewise smooth deformations /„ and a sequence of surface-like subsets of 
the body such that 


g = lim /„ 

n—>-oo 

(1.2) 

G = lim V/„ 

n—>-oo 

(1.3) 

oo oo 


=u n "'p 


n—lp—n 



The limits in ( |1.2| ) and ( |1.3| l are taken in the sense of convergence. A sequence 
n I—^ fn of piecewise smooth, injective functions satisfying ( |1.2| ) and ( |1.3| l is called 
a determining sequence for the pair {g, G) , and each term /„ is interpreted as de¬ 
scribing the body divided into tiny pieces that may individually undergo smooth ge¬ 
ometrical changes and that also may undergo disarrangements, i.e., may separate or 
slide relative to each other. In this context, we write /„ {g, G). From ( |1.3| l we 

see that G captures the effects at the macrolevel of smooth geometrical changes at 
submacroscopic levels, and we call G the deformation without disarrangements. 
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Del Piero and Owen [O proved that, for every structured deformation (k, g, G), for 
every determining sequence n i—/„ for {g,G), and for every point x where g is 
differentiable and where G is continuous, there holds 


lim lim 

r—vO n—^oo 


lj(f„)nB^(x) 


[fn]{y) i^{y) dn (y) 


\Brix)\ 


Vy(a;) - G{x). 


(1.4) 


Here, denotes the (iV—1)-dimensional Hausdorff measure on Br{x) de¬ 

notes the open ball centered at x of radius r, \Br{x)\ denotes its volume (i.e., its 
TV-dimensional Lebesgue measure), J{fn) denotes the jump set of fn, i.e., points 
where /„ can suffer jump-discontinuities, and [/n](j/) ® v{y) is the tensor product of 
the jump [fn] of /„ with the normal u to the jump set. This result permits us to call 
the tensor 

M{x) := Vgix) - G{x) (1.5) 

the deformation due to disarrangements, because it captures, in the limit as n tends 
to infinity, the volume density of separations and slips between pieces of the body as 
described by the approximating deformations /„. We may then regard the tensor 
field M as a tensorial disarrangement density that, for every determining sequence 
n I— fn for (y, G ), reflects the limits of interfacial discontinuities of the approximat¬ 
ing deformations /„. Moreover, ( |1.2| l and ( |1.3| ) along with the definition of M <IID 
yield the alternative formula for the disarrangement density: 


M = \7{ lim fn) - lim V/„. 

n—¥oo n—>-oo 


Consequently, M measures quantitatively the lack of commutativity of the classical 

gradient V and the limit operator lim„_for L°° -convergence. 

The trivial algebraic relation 

Vg = G + M (1.6) 

together with the identification relations ( |1.3D and ( |1.4| l shows that the macroscopic 
deformation gradient V g has an additive decomposition into its part G without dis¬ 
arrangements and its part M due to disarrangements. Because G has invertible 
values, leads immediately to two multiplicative decompositions for V g: 

Vff = G{I + G-^M) = {I + MG-^)G. 


The disarrangement density M and the deformation without disarrangements G 
have an additional property significant in the description of defects and dislocations 
in a continuous body in three dimensions. We consider a smooth surface S with 
smooth bounding closed curve 7 , both contained in a region in the body where g and 
G are smooth. The relation ( |1.6| l and the smoothness of g imply 

0 = ® Wg{x)dx = ® G{x)dx + (p M {x) dx. 

J J -y J -y 

The vector M {x) dx measures the displacement due to disarrangements along 7 
and may be called the Burgers vector p7l for 7 arising from the given structured 
deformation. Application of Stokes’ Theorem to G(x)dx and M (x) dx and use 
of the previous relation yields the formulas for the Burgers vector: 

fM{x)dx= / c\u:\M{x)v{x)dAx = — / cm:\G{x)v{x)dAx. 

J-f Js Js 

The second-order tensor field curlM = —curlG thus determines the Burgers vector 
associated with 7 for every closed curve and corresponds to familiar measures of 
dislocation density (HI 1121 . In this manner, the disarrangement density tensor M 
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determines both the Burgers vector and the dislocation density tensor, both basic 
tools in modelling the effects of submacroscopic defects on the response of solids. 

The tensorial relations ( |1.5| l and ( |1.4| ) yield upon application of the trace operator 
the scalar relation 


lim lim 

r —>-0 n—>-oo 


lj{u)nBr.ix) 


[fn]{y) ■ I'iy) du (y) 


\Brix)\ 


trM{x) 


(1.7) 


in which [fn]{y) ■ >^{y) is the scalar product of the jump and of the normal at y. The 
formula ( |1.7| l tells us that trM is a scalar (bulk) disarrangement density that cap¬ 
tures the components of the jumps of /„ that are normal to the jump set. Moreover, 
this scalar disarrangement density at x, trM{x), allows for cancellation of positive 
and negative contributions of [fn]{y) ■ v{y) at points y near x to the integral on 
the left-hand side of ( |1.7| l. Thus, trM{x) does not distinguish between jumps with 
[fn]iy) ■ v{y) > 0 that pull apart small pieces of the body near x and jumps with 
[fn]iy) • u{y) < 0 that cause small pieces near x to switch places. Because the ap¬ 
proximating fn are injective, the possibility for the case [fn]iy) ■ v{y) < 0 that /„ can 
cause adjacent small pieces of the body to interpenetrate is ruled out. 

Owen and Paroni II13II refined the scalar disarrangement density trM by replacing 
[fn] iy) ■ v{y) by its positive part throughout the jump set of /„, or by its negative part 
throughout the jump set: 


i[fn]{y) ■ J^iy))^ = U\lf^]iy) ■ + IfnKy) ■ v{y)) (1.8) 


[[fn]{y)-v{y)) =\{\[fn]{y)-v{y)\-[fn]{y)-v{y)). (1.9) 

The field ([/„] ■ on the jump set is a scalar (interfacial) disarrangement density 
that measures separations of pieces of the body caused by /„, while the field i[fn]-u)~ 
is a scalar (interfacial) disarrangement density that measures the switching of pieces 
of the body caused by /„. Since we have 

\[fn]{y) ■ u{y)\ = {[fn]{y) ■ Hy))^ + ([/»](?/) • Hy))~, 

the field \[fn]iy) • t^iy)] is a scalar disarrangement density that measures both sep¬ 
arations and switches. We fix a part V of the body, we integrate ( |1.8| l or ( |1.9| l over 
7(/n) n P and use the formula ( |1.7| l to obtain the relations 


liminf [ {[fn]{y) ■ u{y))^dn^ ^{y) 

= i liminf /" \[fn]{y) ■ Hy)\dn^~^{y) 

2 n^oo 

± liminf f [fn]{y) ■ v{y)dn^-^{y) 

2 n^oo 

= i liminf f \[fn]{y)-i'{y)\dn^-^{y)±^ f trM{x)dC^[x). 

2 n-^-oo JJ(^f^)r]V ^ j-p 


( 1 . 10 ) 


Consequently, the limiting behavior of the integral of ([/n](j/) ■ 1^(2/))* in dol l as n 
tends to 00 is determined by the behavior of the integral of \[fn]{y) ■ r'iy)], and we 
restrict our attention to the latter. We expect that 


unlike 


liminf / \[fn]iy) ■ u{y)\dn^ ^(y), 

n->-oo 


liminf [ [fn]iy) ■ >^{y)dn^ ^(y), 
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will depend upon the choice of determining sequence for {g, G). Therefore, we are 
led to consider the most economical manner in which separations and switches can 
arise among the determining sequences for {g, G) : 

:= inf I liminf /" \[f„\{y)-v{y)\(m’^~^{y) :(g, G) I . (1.11) 

The number V' ' (g, G;7^) so defined has the dimension of volume, and we call Vll(g, G; V) 
the (minimal) volume swept out by disarrangements in V for {g,G). If we replace 
I• I everywhere in ( |l.ll| l by ” + ” or everywhere by ” — ”, then we call the number 
V+(g, Gi'P) the (minimal) volume swept out by separations in V for {g,G), and the 
number V~ (g, G; V) the (minimal) volume swept out by switches in V for (g, G). The 
formulas ( |1.10| l imply the simple formulas 

V^{g,G-V) = \v\\g,G-V)±\ [ trM{x)dC^(x) 

Jv 

and, in view of the form of the second term on the right-hand side, raise the follow¬ 
ing basic question: does the volume swept out by disarrangements Vl l(g, G;7^) as 
defined in ( |l.ll| l have an associated disarrangement density which, when integrated 
over V , recovers V' l (g, G; V). If so, what specific information can be obtained about 
the dependence of the integrand upon the structured deformation (g, G)? 

While the setting for structured deformations described in this subsection is quite 
suitable for formulating refined field equations in continuum mechanics [|8j that re¬ 
flect the influence of submacroscopic geometrical changes in a body, this setting has 
not provided answers to questions such as the ones just raised. Part of the difficulty 
with the setting provided in (7!l lies in the choice of smoothness placed on g and its 
approximates /„, while another part lies in the requirement that g and /„ be injec¬ 
tive. An alternative setting provided by Choksi and Fonseca [dl was proposed for 
dealing with such questions and is described briefly in the next subsection. 

1.2. Structured deformations and disarrangement densities in the setting 
of Choksi and Fonseca. We describe here a few essential elements of the treatment 
of structured deformations by Choksi and Fonseca [ Jj. The articles tSjl, 111, m, 
and HSl also provide summaries of that treatment, and HI, |2l, and fl5l provide 
alternative settings for structured deformations. The summary in [3l is intended for 
those interested in immediate applications in continuum mechanics, while HI sets 
the stage for applications of structured deformations to thin bodies IfTIl . The article 
Ifl5l reexamines the results of lH in a broader setting while providing refinements of 
counterparts of the Approximation Theorem and the identification relation ( |1.4| l. 

According to Choksi and Fonseca, a structured deformation is a pair (g, G) in 
which g : D —> , with fl an open subset of the space of N -tuples of real 
numbers, and G : D — R^^^ , with M^^^the space of Ax A matrices with real 
entries. The mapping G is assumed to be integrable on D, G e and 

g is assumed to be in the space SBV{n;R^), i.e., g is a function of bounded varia¬ 
tion with the additional property that its distributional derivative Dg, as a bounded 
measure, has zero Cantor part: 

Dg = \/g£^ + [g]<^un^-\ (1.12) 

Here the integrable mapping Vg is the density of the absolutely continuous part 
Vg£^of Dg with respect to A-dimensional Lebesgue measure , and [g] 0 is 
the density of the singular part [g]®v'H’^~^ Ag with respect to (A—1)-dimensional 
Hausdorff measure . The singular part is concentrated on J(g), the jump set of 
g, and, as usual, [g] denotes the jump in g and v denotes the normal to the jump set 
J(g). It is important to note that Vg in the present setting is no longer the classical 
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gradient of a smooth field and, consequently, need not be curl-free. Nevertheless, 
V g satisfies an integral version of the property of approximation by linear mappings 
that defines the classical gradient of smooth fields . 

Choksi and Fonseca Ul prove a version of the Approximation Theorem with ap¬ 
proximating deformations /„ also in SBV{fl;R^) and with ( |1.2| ) and ( |1.3| l replaced 
respectively by 

/„^ginLi(f!;M'^) (1.13) 

and 


^ fn ^ G weakly in the sense of measures. (1.14) 

We note that no restriction in the form of the accommodation inequality ( |1.1| ) or in the 
form of a requirement of injectivity of g or /„ is imposed in the present context. We 
again use the term determining sequence to describe a sequence n i—/„ satisfying 
( |1.13| l and ( |1.14| l for a given structured deformation [g, G ), and we again write /„ 

{g, G) when ( |1.13| l and ( |1.14| l both hold. The properties of distributional derivatives 
along with relations ( |1.12| l, ( |1.13| l, and ( |1.14| l justify the calculation 

\/gC^ + [g] ® = D lim /„ 

71—>00 

= lim Dfn 

n—>oo 


n—>oo 

= G+ lim {[fn] ^ 


where the convergence indicated in the last three lines is weak convergence in the 
sense of measures. We conclude that the singular parts [/„] 0 of the approx¬ 

imating deformations /„ converge in the same sense and that their limit satisfies 

lim {[fn] 0 = (Vg — G)C^ + [g] 0 (1.15) 


In particular, the restriction of the limiting measure lim„_>.oo([/n] 0 to the 

complement of the jump set J{g) agrees with the corresponding restriction of (Vg — 
G)L^ = M . Consequently, the tensor field M = Vg — G retains in this broader 
setting for structured deformations its identity as a tensor density of disarrange¬ 
ments for (g,G). The formula ( |1.15| l shows that when M = Vg — G ^ 0, while 
all of the measures [fn] 0 are supported on sets J(/n) of -measure zero 

and so have -parts zero, the limit measure lim„_>oo([/n] 0 has £^-part 

non-zero. This observation points to the fact that the jump sets J{fn) can dif¬ 
fuse in the limit throughout the domain fl so that the limiting measure lim„_>oo {[fn]0 
is supported in part on sets of positive -measure. This provides a coun¬ 
terpart in the SBV -setting to the relation ( |1.4| l in which the limit of jumps on the 
left-hand side delivers the -density M. (See HTSl for a detailed derivation of a 
counterpart of ( |1.4| l in a somewhat broader setting for structured deformations than 
SBV.) 

We note briefly that the scalar density of disarrangements trM = tr{Vg — G) that 
counts only normal components of jumps and that emerged in the previous setting 
also appears in the present setting when one takes the trace of every member of 
( |l-15| l: if fn (g, G) , then 

lim {[fn] ■ = tr{Vg - G)C^ + [g] • 


However, as was the case in the setting of Del Piero and Owen, replacement of [/„] • u 
by ([/n] ■v)'^ or by | [/„] • v\ need not yield a limit of the corresponding measures and, 
if a limit exists, the limit may depend upon the choice of determining sequence n i—^ 
fn ■ The setting of Choksi and Fonseca was formulated as a means of resolving these 
difficulties, and we summarize some aspects of that resolution in the next subsection. 
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1.3. Relaxation of energies for structured deformations. Optimal functions 
arising from structured deformations such as the one (fni 

V''(5,G;75) = inf I liminf /" \[fn]{y) ■ i^{y)\dn^~'^{y) :/n'->( 5 ,G)j 

I ” J J(/„)rTP J 

introduced in Section 1.1 can be analyzed using the results of Choksi and Fonseca 
DU on "relaxation of energies" for structured deformations. In that approach, the 
integral \ [fn]iy) ■ t^{y) \ d'H^~^{y) is replaced by an initial energy functional 

Eifn)= [ W{vUy))dC^iy)+ f mn]iy),Hy))dH^-\y) 

JQ. J j{f„)nn 

defined for /„ G SBV{fl;R^). By imposing conditions on the initial bulk energy 
density W and on the initial interfacial energy density tp, the goal is to obtain for 
the relaxed energy I{g, G) defined by 

I{9,G) := inf | liminf ( [ W{S/fn{y))d£^{y)+ 

I \Jn 

[ H[fn]{y),i^{y))dn^~'^{y)) : fn (5,G) 1 

Jj{f„)nn J J 

a representation of the form 

I(g, G)= [ H{yg{y), G{y))dC^{y) + f h{[g]{y), ,.{y))dn^-\y) 

JQ, Jj(g)nQ 

and to deduce properties of the relaxed bulk energy density H and the relaxed inter¬ 
facial energy density h . Because our present interest lies in the case of disarrange¬ 
ment densities, and not on the full energetics of structured deformations, we shall 
restrict our attention to the case W — 0, and we record the following adaptation for 
the case IF = 0 of results from m (see 11131 Theorem 3] for further comments and 
other adaptations). 

Theorem 1.1. Let = {u G : \iy\ = 1}. Let Vibe a bounded open subset of 
and Ip '■ X -G [ 0 , +oo) be such that 

(HI) there exists a constant G > 0 such that 


for all {i,v) G X S^-\ 

(H2) ip{-,v) is positively homogeneous of degree 1: 

tp{ti,v) = ttp{^,u) 


for all t > 0 and (^, u) G x 

(H3) ipf^v) is subadditive, i.e., for all ^ 1,^2 S and v G , 

^>(6 + 6 , V) ^ tp{^l,u) + tp{^2,u). 

Then, for any p > 1, if we define 

/(g, G) ;= inf I liminf [ 'ip{[un],v) d'H^~^ : Un G 

L Jj(Ur,)nQ 

Un ^ g in L^{n; K^), ^ G, 

sup (|VM„|ip(o.]RNxJV) -f \DUn\{H)) < -too L 
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we have 


I{9,G) 


where 


Hiyg(x), G{x)) dC 


N 


' J(s)nn 


hi[9Kx),v{x))dH^ ^{x 


H{A, B) ■- inf 


'ip{[u],v) dH 


N- 


^ :uG SBV{Q-R^), 


j j{u)r\Q 

u\dQ = Ax, |Vu| e LP{Q), 


VudC^ 



(1.16) 


and 


:=inf ( f ^Ij{[u],v) dU^-^ : u G SBViQr,;R^), 

I J j{u)r\Qrj 

■“IsQ, = Vu = 0 a.e.|, 


(1.17) 


with 


U^,r,{x) 


0 if — ^ ^ X ■ r] < 0, 
^ if 0 ^ X ■ r] < ^. 


(1.18) 


Here, Q = (—1/2,1/2)^ and Qr^ denotes the unit cube centered at the origin and with 
two faces normal to rj. 


In the right-hand side of ( |1.17| l we have corrected an inconsequential misprint that 
is present in the corresponding formula in Theorem 3 of II13I . 

Another approach to relaxation of energies for structured deformations in the full 
BV setting is provided in DTJ. A structured deformation in [U is a pair {g, G) G 
BV'^{Vt-W^)xBV{n-R^^^),where BV'^{H,W^) := {u G BV{n;R^) : Vu G BV{n;R^^^)}. 
The counterpart of the Approximation Theorem in this context asserts that there ex¬ 
ists a sequence /„ G BV'^{H,R^) such that both f„ ^ g and V fn -g G in the 
H -norm. In this case we write fn-^{g,G). 

The energy functional considered in [TJ, under assumptions on the initial bulk and 
surface energy densities similar to the ones in 14], reads 


E{fn) = f WiyUy),y^U{y))d£^y+ [ mn]{y),Hy)) 

Jn JJUA 

+ / M[^fn]iy),Hy))dn^-\y), 

and the relaxed energy I{g, G) is defined by 

I{g,G) :=inf|liminf£;(/„) (g,G)|. (1.19) 

A crucial result in [11 is that ( |1.19| l can be divided into two first-order relaxed 
energies, namely, I{g,G) = Ii{g,G) + hiG), where the term Ii{g,G) captures the 
structured deformation, whereas / 2 (G) only depends on the deformation without 
disarrangements G. In the relevant case for the present paper, i.e., W = ipi = 0, the 
results in [Ti| give I 2 = 0 and 


Ii{ 9 ,G) := inf 


liminf [ 'tlj{[fr,]{y),u{y))dn^ ^{y) ■ fn {g,G) 


Defining SBV^{n;R^) := {u G SBV{n-R^) : Vu G 5'W(D;, the following 
representation theorem holds 






FORMULAS FOR RELAXED DISARRANGEMENT DENSITIES 


9 


Theorem 1.2 (see Ifll Theorem3.21). Forevery {g,G) G SBV‘^{n-,R^)xSBV{Q\R^^^), 
given ij} under the same hypotheses (H1)-(H3) ofTheorem \l.l\ we have that 


I{9,G)= f H{G{x) -Wg{x))dC^ + f h{[g]{x),n{x)) dH^ ^(x), 

JQ J J{g)r\Q, 

iven A G ^ G and g G , 

H{A):=mi\ [ ^{[u],v) dB^-^ : u € SBV'^iQ;R^), 

I J j(u)nQ 

u\dQ = 0, Vm = A a.e. in 

(i(^, 77 ) := inf I f ip{[u\,u) d'H^~^ : u e SBV'^{Qrf',R^), 

I J j(u)nQr) 


' j{u)r\Q 

— 0 Ct.6. IIT Q ^ .J 

with u^^r] defined as in CTsl i. 


( 1 . 20 ) 


( 1 . 21 ) 


Remark 1.3. It is worth noticing that the minimum problems defining ( |1.20| l and 
( |1.21| l are formally performed in SBV'^{fl; R^), but the result is the same if SBV"^ is 
replaced in these relations by SBV, due to the requirement that Vu be constant. 


1.4. Explicit formulas for relaxed disarrangement densities. Owen and Pa- 
roni |[l3l applied Theorem |l.l| to the specific disarrangement densities \[fn]iy) ■ v{y)\ 
and ([/„](?/) • v{y))^ introduced in Section 1.1 and obtained for each of these densi¬ 
ties an explicit formula for the corresponding relaxed disarrangement densities H 
in ( |1.16| l and h in ( |1.17| l. Among their results ( flSlI .Theorem 4) are the following 
(obtained by setting L{x) = / in their Theorem 4): 

Theorem 1.4. The initial disarrangement densities 

( 1 - 22 ) 


satisfy the hypotheses (H1)-(H3) in Theorem 1.1 and have relaxed disarrangement 
densities given by 


H^-^{A,B) = \triA-B)\, = (1.23) 

and 

H^{A,B) = itriA-B))^, h±(C,u) = 


Specifically, when the minimal volume swept out by disarrangements V' ' (y, G; V) 
is defined in the Choksi-Fonseca setting by ( |l.ll| l, then ( |1.23| l yields the explicit for¬ 
mula 

P''( 5 ,G;P)= [ ItriVgix) - G{x))\dC^{x)+ [ \[g]{x) ■ v{x)\dH^~'^{x) (1.24) 

Jv Jj{g)nv 

for the (minimal) volume swept out by separations and switches among approxima¬ 
tions fn that determine {g, G). Relation ( |1.24| l provides answers in the setting of 
Choksi and Fonseca to the questions raised at the end of Section 1.1: pl l {g, G; V) has 
both a bulk disarrangement density \tr{Vg — G)| = |frM| and an interfacial disar¬ 
rangement density [[g] • v\. Similarly, Theorem |1.4| shows that the (minimal) volume 
swept out by separations alone, V+(g, G;P), has the bulk disarrangement density 
(trM)+ and the interfacial disarrangement density {[g] • , with a corresponding 
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result for V~ {g,G-,V), the (minimal) volume swept out by switches and interpene¬ 
trations (the approximations /„ in the Choksi-Fonseca setting are not required to 
be injective, so that interpenetrations can arise there, unlike in the setting of Del 
Piero-Owen). 


1.5. Summary of the research presented in the present article. In the proof 
of Theorem 1.4 given in II13II . the significant part of the argument addresses the 
verification of the inequality 


if'•1(^,5) ^ \triA-B)\ 


(1.25) 


where ifl l (^1, B) is given by the right-hand side of ( |1.16| l with ^p{[u],l^u) replaced by 
= |[m] • z/|. This inequality was proved in IIT31I by constructing a family 
Us of piecewise affine mappings on the unit cube Q each of whose jump set J{us) 
is formed by two (planar) ends and by a lateral surface constructed from solution 
curves of the differential equation x = {A—B)x. The lateral surface, by construction, 
contributes nothing to the integral |[u] ■ v\ , and the contributions of 

the two ends can be calculated explicitly for A — B lying in a dense subset of . 

Proposition 5.2 of 01 provides sufficient regularity of H^'\A,B) to establish ( |1.25| l 

for all A-BgR^^^. 

As one of the main results in this article, we provide an alternate, shorter proof 
of ( |1.25| l that employs a different family Us of piecewise affine mappings that does 
not involve solution curves of i: = (A — B)x. Our approach is based on the following 
observation. With A,B G , p> 1, and with Q = (—1/2,1 /2)^ there hold 


|tr(A — i3)| < inf I f \[u\{x) ■ iy{x)\: u G SBV{Q;R^), 

I Jj(u) 

u{x) = Ax on dQ, Vu G L^{Q), J Wu{x)dC^(x) = sj. 
^ inf I f \[u\{x) ■ i^{x)\d'H^~^{x) : uGSBV{Q;R^), 

L J Jiu) 


(1.26) 


I J{u) 

u{x) = 0 on dQ, Vu = B — A a.e. ^. 


The first follows by moving the absolute value outside the integral and using the 
Gauss-Green Theorem for the space SBV{Q;R^) of special functions of bounded 
variation, while the second follows by noting that if u satisfies the last set of condi¬ 
tions, then the function x i—^ u{x) + Ax satisfies the first set of conditions. In this 
paper, we wish to show that 


inf 


\[u]{x)-v{x)\dH^-^{x) : uG SBV{Q-,R") 


N\ 


' J{ii 


u{x) = 0 on dQ, Vu = B — A a.e. > < \tr{A — B)\ 


(1.27) 


so that the two infima in ( |1.26| l have common value \tr{A-B)\. 

The second main contribution of the present research concerns the alternative 
approach to structured deformations and to relaxed energies due to Baia, Matias, 
and Santos [T] discussed at the end of Subsection [O According to that discussion 
the second infimum in ( |1.26| l (see ( |1.20| l and Remark |1.3| l 


inf 


' J(h 


\[u]{x)-v{x)\dH‘^-\x) : u G SBV{Q; R‘^) 


u(x) = 0 on dQ, Vu = B — A a.e. in Q 
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is the bulk disarrangement density for the same interfacial disarrangement density 
^11 ([u], ^1.22) studied by Owen and Paroni in the setting of Choksi and Fonseca. 

Consequently, our proof of ( |1.27| l establishes the equality of the bulk disarrange¬ 
ment densities obtained in two different settings for structured deformations. Thus, 
the geometrical significance of the expression \tr{A — B)\ described in I113I . namely, 
a volume density of volume swept out by non-smooth, submacroscopic geometrical 
changes, is strengthened by the fact that one and the same expression arises from 
two different schemes of relaxation. We note that the two different schemes of re¬ 
laxation also deliver the same formula for the (relaxed) interfacial disarrangement 
density h: h = -011 (see 1131 for the routine verification that applies to both schemes). 

The explicit formulas for disarrangement densities considered here in the context 
of structured deformations will provide scalar fields that can enter as variables in 
constitutive relations for the response of three-dimensional bodies. For this purpose, 
frame-indifferent variants of the specific fields obtained here are available through 
known factorizations of structured deformations in which the factor that tracks dis¬ 
arrangements is unchanged under changes in frame [7]. Our explicit formulas also 
are starting points for the study of examples in other contexts involving structured 
deformations: second-order structured deformations II 1411 in which second gradients 
and their limits enter into submacroscopic changes in geometry, as well as processes 
for dimension reduction I 111 in the presence of disarrangements that describe thin 
structures undergoing submacroscopic slips, separations, and switches. 

In Section 2 we provide a "tilted cube" construction for the family of functions 
employed in proving ( |1.27| l. The common orientation of the tilted cubes is determined 
in Section 3 by means of a known result on the isotropic vectors of symmetric linear 
mappings. The proof of ( |1.27| l is completed in Section 4, and the paper concludes in 
Section 5 with some additional explicit formulas for disarrangement densities. 

During the review of this article, the research 1161 was brought to our attention in 
which explicit formulas for the bulk and interfacial relaxed energies are established 
for a broad class of purely interfacial initial energies that includes the ones studied 
here. 


2. Proof of the upper bound inequality 

In what follows, a proof of ( |1.27| l is given. The proof requires the following instance 
of Lemma 4.3 in IflOl . 

Lemma 2.1. Let M e and a bounded open set Q C be given, with Q 

having Lipshitz boundary. There exist a number C{N) > 0, independent of M and 
Q, and u G SBV{n]M.^) such that 

(1) u|an = 0 

(2) Vu = M, — a.e. on Q 

(3) |D"u|(Q) ^ C(iV)||M||/:^(Q). 

Here, Vu and D^u denote the absolutely continuous and the singular parts of 
the distributional derivative Du = VuC^ + D^u of u, and \D‘u\ denotes the total 
variation of the singular part. In addition, ||M|| := is the Euclidean 

norm of the matrix M. We shall now use the Lemma to verify ( |1.27| l for M = A — B. 
To this end, let an integer n ^ 1 be given and consider the frame 

d^n-.= Q\{l-^)Q. 

We may apply the Lemma to obtain an SBV function : Tn satisfying 

• =0 

• =M, - a.e. on Tn 
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• the total variation | of satisfies 

[ |[«(-)]|(a:)dH^-i(x)<C(7V)||M|| (i-( 1 - 45 )^) ( 2 . 1 ) 

In preparation for defining an appropriate function u on Q\Tn = (1 — we 

write M := \{M + M^) for the symmetric part of M, and we choose an orthonormal 
basis Ci, i = 1,..., iV of that consists of eigenvectors of M: 

Mci = AiCi, i = 1,... ,iV. 

We let m be a positive integer and cover ( 1 — ;^)Q by a collection Cn^m of congruent, 
non-overlapping open cubes k = each of edge-length 1 /m, each 

with the pair of opposite faces orthogonal to the unit vector Rci , for i = 1,..., iV. 
Here, R is an orthogonal NxN matrix, RR^ = R^R = / , to be determined presently. 
We require in addition that each each cube satisfies 


We denote by ^ the center of , and we define Un^m '■ (1 — 

„ /..'i ._ifa;G ( 1 -;^)QnC'*„ 

\o if a; G (1 - ;^)Q\ 


for some fc = 1 ,..., Kn,m, 

/^k 

^n,m’ 


Using standard reasoning we conclude that Un,m G SBV{{1 — with 

= M, — a.e. on (1 — :^)Q. Moreover, the trace of Un^m on d{{l — ;^)Q) 

is bounded pointwise by ^ ||M||. Consequently, the function u'm '■ Q defined 

by 




(a;) := 


fora;GJ-'„, 


for X G (1 n-t-2^^ 

belongs to S'ilU((5; K^), has gradient M, C^— a.e., and has zero trace on OQ. More¬ 
over, the jump set of Um^ satisfies 


C U 9((1 - ^)Q) U J{Un,m). 

Since Um'’ has outer trace 0 on 9(1 — :^)Q, there holds for 


( 2 . 2 ) 


— a.e. X in 


and, consequently, 


()](x)|<^||M|| 


y/N 


We note from ( |2.1| l that 


[Wm^](a;) ■ v{.x) dn^ i(x) < ||M|| 2N{1 - 


2 \W-1 


/ 

• v{x) 




and it remains to obtain a corresponding estimate for ^ 

To this end, we note that 


n+2J h 
[U'm\{x) ■ v{x) 


(2.3) 

(2.4) 

dU^-^ix 


K„ 


J{Un,m) C U dCl^ , 
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and we shall seek an upper bound for [um'^]{x) ■ v{x) ^{x). For 

each fc = 1 ,..., Kn,m and i = 1,..., we denote by and (j^n^m the two faces 
of the cube C^ „j S Cn,m orthogonal to Rci. We note that one face 0^’*j+of has 
outer normal u+j, = +R , while the opposite face (t)n'm has outer normal 


R 6 ?;. 


We suppose now that the face ^n,m ^ satisfies 

C (1 - ^^)Q. 

yk' 


'^i,k 


(2.5) 


Then there is a cube ^ e C„ „i that shares the given face with ^, and we have 
at each point x G 


Afc.*+ 


= iM{x-ct,m)-M{x-c^ ,^))-iy+^{x) 


— Cn,m) ' 

= - MRci ■ Rci 

m 


so that 


■ vix) 


dn 


N-1 


(x) = 


1 


j,k,i+ m 


\MRei ■ Rei\d'H^ ^(a:) 


1 


( 2 . 6 ) 


= — j^\MRei ■ Rci 


The same argument shows that if 


ik.i— 


’.m C (1 - 


then 


[ [u^^'>]{x) ■ v(x) dU^ ^{x) = ^jrlMRei ■ Re^ 


(2.7) 


If ( |2.5| l holds for i = 1,..., iV, then we may sum the last relation over i to conclude 
that 

I [ul^^Kx) ■ I^{x) dH^ i(x) IMiiei • iicil 


i=l 

N 


>- 


iN 


iN 


MRci ■ Rt 


i=l 

N 


( 2 . 8 ) 


YR^mr 


Ci * Ci 




= ^MR^MR)\ = ^\trM\. 

In ( |2.8| ) equality holds if and only if all of the numbers MRci ■ Rci , i = l,...,iV, 
have the same sign: 


[MRci ■ Rci ){MRej ■ Rcj) ^0, for i, j = 1,..., TV. 

The last two inequalities lead us to consider the problem 


(2.9) 


N 


Find min \MRei ■ Rei\ ^ |trM| = |trM|, 

RHY—I . 

1—1 


( 2 . 10 ) 
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with equality holding if and only if there exists an orthogonal matrix R satisfying 


3. Aside on isotropic vectors 


We note that the sign inequality ( |2.9| l suggests looking for unit vectors v such that 

Mv ■ V = 0, (3.1) 

the isotropic vectors for M ||5l. In particular, in the special case tr M = 0, the 
existence of N mutually orthogonal isotropic vectors vi,... ,vn would insure that the 
matrix R defined hy Rci = Vi for i = 1,... ,N would satisfy ( |2.10| l in the form 0 = 0. 
More generally, even when tr M ^ 0, the existence of isotropic vectors is useful. In 
fact, the symmetric matrix M — ^{tr M)I has zero trace, so we suppose that there 
exist N mutually orthogonal isotropic unit vectors vi,... ,vn (or M — j^{tr M)I . The 
relation ( |3.1| l with M replaced hy M — ^{tr M)I then becomes 

0 = - ^{trM)I^Vi ■ Vi = Mvi ■ Vi - 


so that Mvi ■ Vi = for i = 1,..., A. Again, if we define a linear mapping R on 
K^by Rci = Vi for i = 1,..., A then R is orthogonal, it satisfies the sign inequality 


for M 


, and it delivers equality in ( |2.10| l in the form X)i=i 


tr M 


N 


= \tr M\. 


The following result (|hi|. Corollary 15) provides the desired existence of complete 
orthonormal sets of isotropic vectors. 


Theorem 3.1. A symmetric matrix A e R^^^ possesses an orthonormal set of A 
isotropic vectors if and only if trA = 0. 


This theorem and the preceding discussion permit us to conclude: for every matrix 

M e 

N N 

min iMRci ■ ReA = min IMRci ■ ReA 

RRr=I^' 

1—1 1—1 

= \tr M\ = \trM\ , 

and a minimizing rotation matrix R is one carrying the orthonormal basis of R^ con¬ 
sisting of eigenvectors of M into an orthonormal basis of R^ consisting of isotropic 
vectors of M — T(<r M)/. For this minimizing rotation matrix, we have 

\MRei ■ ReA = ^ \trM\ for i = 1,..., N. (3.2) 

We remark that minimizers are not unique, in general, even when one eliminates 
trivial permutations of isotropic vectors. In fact, for A = 3 there are examples of 
minimizers for which two of the three terms in \MRci ■ ReA vanish, while the 
third equals \trM\ , so that only two of the three vectors Ret are isotropic vectors for 
M . 

For the convenience of the reader, we provide the recursive step used in proving the 
existence of orthonormal bases made up of isotropic vectors for a traceless symmetric 
matrix A e R^^^. We interpret A in the usual way as a linear mapping on R^, 
endowed with the standard inner product. Then the nullspace KerA of A and 
its orthogonal complement {KerA)^ are complementary A-invariant subspaces of 
R^, and all vectors in KerA are isotropic vectors for A. If {KerA)^ is the zero 
subspace, then A = 0 and every vector in R^ is an isotropic vector for A, and every 
orthonormal basis of R^ meets the desired requirement. If {KerA)-'- is not the 
zero subspace, then we seek additional isotropic vectors for A in (KerA)-'- . To this 
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end, the traceless symmetric linear mapping A ^ 0 has both positive and negative 
eigenvalues so that 


min Au • u < 0 < max Au • u 

|u| = l \u\ = l 


and, since the unit sphere in is connected and since the quadratic form u i—s- 
Au ■ u is continuous, there exists a unit vector vi € such that Avi • vi = 0. Writ¬ 
ing vi as a sum of two orthogonal vectors, one in KerA and the other in (KerA)-^ 
and using the invariance of (KerA)-^ under A shows that we may without loss of 
generality assume that vi e (KerA)^ . The linear span Lsp{KerA U {vi}) has di¬ 
mension one larger than that of KerA and consists solely of isotropic vectors for A. 
Consequently, we need to search for isotropic vectors of A in {Lsp{KerA U {I’l}))-'- 
which has dimension one less than {KerA)^ . To procede further, we define a linear 
mapping Ai on by 


Ai = A — vi ® Avi — Avi 0 vi 


where the formula (a ® b)v := {b ■ v)a , for all a,b,v G , defines the standard tensor 
product a 0 5 G Lin{R^;R^) . From the fact that vi is an isotropic vector for A 
and from the formula tr((a 0b)=a-bitis easy to see that Ai is traceless; because 
(o(g)6)^ = big)a, it follows that Ai is symmetric. In addition, if v e {Lsp{KerAu{vi}))-^ 
is an isotropic vector for Ai, then we have not only v • = 0 but also 

0 = Aiv ■ V 

= {Av — {Avi • v)v\ — (ill • v)Avi) • V 
= Av ■ V — {Avi • v)(vi ■ v) — (vi • v)(Avi • v) 

= Av ■ V. 


Thus, every isotropic vector for Ai that is in (Lsp{KerA U {vi}))-'- is an isotropic 
vector for A, and (i\m{{Lsp{KerA\J{vi}))^) = dim((iferA)-*-) — 1. To be able to apply 
the forgoing considerations to Ai, we need only show that {Lsp{KerA U {I’l}))''' is 
invariant under Ai. To this end, let v G {Lsp{KerA U {vi}))-'-, G KerA, and 
a G M be given, and consider 


Aiv ■ (Vf, + avi) 


Aiv ■ Vk, + Aiv ■ avi 
V • AiVk + cxv • Aivi 

0 -I- ar • {Avi — {Avi 0 vi)vi — (vi 0 Av\)vi) 
av ■ {Avi — (fi • vi)Avi — {Avi • 
av • (Avi — Avi — 0) = 0. 


We may conclude that Aiv G {Lsp{KerA U as desired. In the third line of the 

above computation we have used the side-calculation 


V ■ AiVft = V ■ {A — vi ^ Avi — Avi ® vi)vk, 

= V ■ Av,^ - {Avi ■ ■u„)(r • wi) - {vi ■ w„)(z; • Avi) = 0. 


The first term on the last line vanishes because v^ G KerA, the second vanishes 
because v G {Lsp{KerA U {vi}))-'- and the third vanishes because vi G (KerA)^. 
The search for isotropic vectors for A on the A-invariant subspace (KerA)-^ may now 
be replaced by the search for isotropic vectors for Ai on the Ai -invariant subspace 
{Lsp{KerA U {wi}))^ of dimension one less than that of {KerA)-^. 
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4. Completion of the proof of the upper bound inequality 
We may use ( |3.2| l and the formulas ( |2.6| l, ( |2.7| l to conclude that: if has a face 


f^,i± 

'^n,7n 


^ (1 - T^)Q’ then 


■ v{x) 


dn 


\trM\ 


Nm^ 


N 


(4.1) 


C (1- iv^)Q, 


On the other hand, if a face (l)n!m of C'n.m ^ Cn,m fails to satisfy 
then the argument used to verify ( |4-l| l may be applied to '/'Sn (1 - ^)Q to 
conclude that 






dn^-\x) < 


(4.2) 


We now consider the cube e Cn,m and choose , one of its 2 ^ vertices. 
Exactly N faces ,j = I,..., N, of meet at Because each cube C* ^ G 

Cn^m for k = 1,..., Kn^m cuu be obtained from ^ by a unique translation , the 
choices „j and ^ induce via Tk an assignment of N faces , j = 1, ..., iV to 
■ It is easy to show that for all k,k' = 1,..., Kn,m 

k'^k^ : j = 1 ,.. . ,iV} n : J = 1,... ,7V} = 0, 

i.e., the set of N faces assigned to different cubes are disjoint. If we now apply the 
mapping 


to each cube in the collection 

C int J 

m ’TYt ' S I—^ 1 




k 

n.m 


a r ‘ n 

^ ^n.m • 


k 

n.m 




then all of the faces so obtained will be included in (1 — j^)Q, and we may 
apply ( |4-1| > to each such face to obtain for each e 

N 

§4- 


• iy{x) 


dn 


N-1 


(x) = = \trM\C^{C’^^J. 


We may sum both sides over the cubes G to obtain 


N 


Ck gCint j = ’’ 




dn 


N-lt 


{x) = \trM\C^{\Jck_^c^..Cl^). 


The faces represented on the left hand side need not include all of J{un,m) C ^ , 

because some faces of cubes G that are also faces of cubes G 
are left out, while proper subsets (1 ~ faces ^Iso are left out. 

However, for those parts of J{un,m), we may use ( |4.1| l and ( |4-2| ) to estimate the in- 

dn^~^{x ), along with the fact that the cubes 


tegrals / 

■ v{x) 

hi':rPr\(i-:^)Q 



whose faces contain these parts of J{un,m) all must contain points of 9(1 — 
and must together cover 9(1 — ;^)Q. Combining all of these contributions to 




u!pp'^]{x) ■ ^{x) dn^ ^(x) we obtain 



■ v{x) 




,6Ci' 


Ct 


^2 |trM| £^(Uc._„6(c„.„\c-„)C((,,„)- 


(4.3) 
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The factor of 2 = ^ in the last expression reflects the fact that the -measure of 
some of the cubes in the collection Cn,m\C™%i has been counted more than once but 
no more than 2N times through the use of the bound ^A.2) . The relations ( |4-3| l, ( |2.2| ), 
( |2.3| l, and ( |2.4| l now yield the relation 




y/N 


+ C{N)\\MUl-il-^f). 

We use in turn ( |4.4| l to obtain an upper bound for 

[ ,, {u’m^](.x)-v{x) dn^-\x). 


2 \N-1 


(4.4) 


Let e > 0 be given and choose n so large that C{N) \\M\\ (1 — (1—;^)^) < £ and, for 
such an n, choose TO so large that ^ ||M|| 2 iV(l—Because (1 — 
has finite measure, we may choose to larger if necessary so that the cover Cn,m 
of (1 “ satisfies ^GC„,„C'n.m) < >C^((1 - :^)Q) + e < 1 + e. Finaliy, 

because 9(1 — :^)Q has zero measure and is covered by we may 

again choose to larger, if necessary, so that 2 |trM| {VjQk ^g(c„ Cim) < e. 

We conclude that for n and to so chosen 


/ .. 





and, since £ > 0 was arbitrary, that ( |1.27| ) holds. □ 


5. Additional explicit formulas for disarrangement densities 

Our discussion above shows that the particular choice of interfacial measure of 
disarrangements 

/ \[u] ■ v\ dH^ (5.1) 

J j {u)r\Q. 

for deformations m of a region 0 c leads in both the Choksi-Fonseca relaxation 
scheme UJ and in the Baia-Matias-Santos relaxation scheme |T1| to one and the same 
bulk density of disarrangements 

[ \tr{\7g-G)\dC^ 

Jn 

for structured deformations {g, G) of that region. Moreover, our analysis here pro¬ 
vides an alternative to the proof of this result given in II13I . In that article, it 
was observed that replacement of |[m] ■ v\ by its positive part {[u]- v)^ = KlM ' ^1 + 
[u] • v) results in the replacement of \tr{Vg — G)] by its positive part {tr{Vg — G))+ = 
l(|tr(V 5 ~ G)| + tr{Vg — G)) in the relaxed bulk disarrangement density. (An anal¬ 
ogous result holds for the negative parts, obtained by replacing " +" by " —" in the 
definition of the positive parts.) As pointed out in IIT31 . {tr{Wg — G)(a;))+ may now 
be interpreted as the minimum volume fraction at a point x € fl that can be swept 
out by submacroscopic separations associated with deformations approximating 
the structured deformation {g, G). Moreover, {tr{Vg — G){x))~ is the minimum vol¬ 
ume fraction at x swept out by submacroscopic switches and interpenetrations, so 
that \tr(yg — G)(a:)| = {tr{Vg — G)(a;))+ -t {tr{Vg — G){x))~ is the minimum volume 
fraction swept out by submacroscopic separations, switches, and interpenetrations. 










18 


ANA CRISTINA BARROSO, JOSE MATIAS, MARCO MORANDOTTI, AND DAVID R. OWEN 


The presence of the inner-product [u] • v in the initial interfacial density ( |5.1| l tells 
us that only normal components of jumps will contribute and that alternative ini¬ 
tial interfacial densities are required in order to capture contributions of tangential 
components of jumps. In the remainder of this section we shall provide alternative 
initial interfacial densities that not only capture contributions of tangential compo¬ 
nents of jumps but also lead to specific formulas for the relaxed bulk disarrangement 
density via the "tilted cube" construction provided in Sections 2 and 4 above. 

Let a € be given and consider the following replacement for ( |5.1| l 

[ \[u]-al dn^-^ (5.2) 

J j {u)no> 

in which the normal component [m] • v of the jump in u is replaced by the component 
[u] • a in the direction of a. To follow again the relaxation scheme in [iTl we let A,Bg 
^NxN given and require not only u G SBV[Q, but also 

u\aQ = 0 ) Vm = B — a, — a.e. in Q. (5.3) 

We now may use the Gauss-Green formula and ( |5.3| l to write 


a\dn^-^ = 

f 

l([« 

■ a])v\ dR^~ 

1 



lj{v 

)nQ 




i 

([u 

a])vdR’^~^ 




Jj{ 

u)r\Q 




= 


f V{u- 

a)dC^ + [ 

{u 

a)vdR^ ^ 



Q 

JdQ 


= 


f {wuy 

adC’^ + / 

(0- 

a)vdR^~^ 



Q 

JdQ 


= 

{B 

— Ay a 





For the "tilted-cube" construction provided in Sections 2 and 4, we replace the matrix 
M hy B — A , and the relation ( |2.6| ) has here the following counterpart 

[ [ul^^]{x) ■ a dn^-\x) = [ ([iiM](x) ■ aMx) dn^-\x) 

= [ -\{{B- A)Re^ ■ a)ReA dR^-\x) 

= —V \[Rei- {B - A)'^a)Rei \ , 

^iV I I 

and this formula leads to the following counterpart of ( | 2 . 8 | l: 

N N 

51 /, [ut\x)-adR^ ^{x)=—^\{Rei-{B-AY'a)Rei\ 


2=1 

N 


>- 


m 

1 


N 


'^^{Rci ■ {B — A)^a) Rci 


i=l 


(5.4) 


,7V 


\{B-Afa\ . 


The method employed in Sections 2 and 4 (where the symbol M was used in place 
of B — A) then requires the choice of a rotation R for which equality holds in the 
second line of ( |5.4| l. If {B — Ay a 7 ^ 0 we may choose R to be any rotation satisfying 
Rei = {B — Aya/ | (B — A)^a|, and this requirement is then met, because {Rei-(B — 
Ay a) Rci = 0 for i = 2,..., N. It {B — Ay a = 0, then R can be chosen arbitrarily, 
for example, R = I suffices. 
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These observations show that the analysis in Section 4 for ( |5.1| l may be carried 
out step by step for the alternative initial density ( |5.2| l, provided that we replace 
everywhere in Section 4 \trM\ = \tr{B — A)| by \{B — A)'^a\, the Euclidean norm of 
the vector {B — A)^a. If we now define 

H{A,B,a)--hA\ f \[u]ix) ■ a\dn^-\x) : u € SBViQ-,R^), 

I Jj{u) 

u |aQ= 0, Vm = B — a a.e.j’, 

then our observations amount to the formula 

H{A,B,a) = \{B- Afa\ (5.5) 

for the relaxed bulk energy density corresponding to the initial interfacial energy 
( |5.2| l and arising from the scheme [T]. Moreover, an argument similar to that used in 
establishing ( |1.26| l shows that the formula ( |5.5| l also holds for the relaxed bulk disar¬ 
rangement density according to [4]. In the context of a given structured deformation 
(ff, G) on a region ft, ( |5.5| l implies that the particular choice of initial interfacial dis¬ 
arrangement 

/ \[u] ■ a\d'H^~^ 

J j {u)nfi 

for deformations m of a region fl c leads in both the Choksi-Fonseca relaxation 
scheme Ul and in the Baia-Matias-Santos relaxation scheme fl’l to one and the same 
relaxed bulk disarrangement density 


[ \{Vg-Gfa\dC^ 

Jn 


(5.6) 


for structured deformations {g, G) of that region. The integral in ( |5.6| l represents the 
most economical way of introducing jumps in the direction of a while approaching in 
the limit the given structured deformation {g, G), including both jumps normal and 
tangential to the discontinuity surfaces of approximating deformations u. 

We note also the formula 


max H(A,B,Si) = \\B — A\\ 


row max 


where on the left di,...,dN denotes the standard basis of and on the right 
\\B — ^11 row max deuotes the maximum of the Euclidean norms of the rows of B — A. 
The mapping || jlrowmax ■ —S' K turns out to be a norm on , and our 

interpretation of the integral in ( |5.6| ) leads us to interpret the integral 


||(V5-G)(x) 


I row max 


dC^ (a;) 


row max 


as a bulk measure of disarrangements that takes into account at each x € Q the di¬ 
rection that maximizes the relaxed bulk energy densities H{Vg{x), G{x), 6i) for 
i = I,..., N. The bulk disarrangement density maxi=i^,.._Ar if(A, B, S^) = \\B — A\\ 
satisfies 

max H{A, B, Si) ^ ini < max [ \[u]{x) ■ Si\d'H^~^{x) : 

z=l....,7V ' 

ue SBV{Q;R^),u\aQ = 0, Vm = S-Aa.e.j, 

and need not be the relaxed bulk energy density corresponding to the initial interfa¬ 
cial energy max / \[u]{x) ■ Si\dH^~^{x) . 
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